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Abstract 

Usually one finds that dissipation tends to make a quantum system more classical in nature. In 
this paper we study the effect of momentum dissipation on a quantum system. The momentum of 
the particle is coupled bilinearly to the momenta of a harmonic oscillator heat bath. For a harmonic 
oscillator system we find that the position and momentum variances for momentum coupling are 
respectively identical to momentum and position variances for spatial friction. This implies that 
momentum coupling leads to an increase in the fluctuations in position as the temperature is 
lowered, exactly the opposite of the classical like localization of the oscillator, found with spatial 
friction. For a parabolic barrier, momentum coupling causes an increase in the unstable normal 
mode barrier frequency, as compared to the lowering of the barrier frequency in the presence of 
purely spatial coupling. This increase in the frequency leads to an enhancement of the thermal 
tunneling flux, which below the crossover temperature becomes exponentially large. The crossover 
temperature between tunneling and thermal activation increases with momentum friction so that 
quantum effects in the escape are relevant at higher temperatures. 

PACS numbers: 



I. INTRODUCTION 

The effect of spatial dissipation on the classical |l| and quantum ^ dynamics of a system 
is well understood. On a microscopic level, dissipation arises from bilinear coupling of the 
system coordinate to the displacement coordinates of a harmonic bath. Classically the bath 
modes obey forced oscillator equations of motion, which may be solved formally in terms of 
the motion of the system. These are then inserted in the system equation of motion which 
then takes the form of a generalized Langevin equation. 

n 

Over twenty years ago, Caldeira and Leggett [3] took advantage of this equivalence to 
study the effect of a dissipative bath on the quantum dynamics of the system, paying special 
attention to the quantum tunneling effect. Their central conclusion was that dissipation 
reduces the tunneling probability, however, it does not destroy it completely. Hence, the 
possibility of observing macroscopic quantum tunneling. 

The detrimental effect of an interaction between the system and its environment on 
quantum phenomena makes intuitive sense. Consider first the localization of a particle in 
space. It is well known [4] that the position variance of a dissipative harmonic oscillator 
becomes smaller as the dissipation strength is increased. In the limit of very large Ohmic 
friction, the bath can localize the particle completely, without violating the uncertainty 
principle |5)]. The bath may be thought of as creating an effective particle with a very large 
mass, and such a heavy particle may be localized. 

One has a similar picture of how the environment destroys tunneling. For a dissipative 
parabolic barrier, it is well known that diagonalization of the system-bath Hamiltonian leads 
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to an unstable mode, whose frequency decreases as the coupling strength increases [a, 17|. 
Tunneling occurs by transmission through this collective unstable mode. Since its frequency 
is smaller, it is a broader barrier, the action needed to cross it increases and the tunneling 



probability decreases 



The same qualitative picture holds at low temperatures below 



the crossover temperature separating the tunneling and activated barrier crossing regimes 
[2|, y, l9[. Spatial dissipation also reduces the crossover temperature [I9j, it is proportional to 
the collective mode barrier frequency. This lowering also fits in with the general observation 
that dissipation causes quantum systems to behave more like classical systems 

All of these conclusions are based on the extensive study of the quantum dynamics of 
dissipative systems, where the Hamiltonian can be brought to the form of bilinear coupling 
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between the system and bath coordinates. There is a quahtative difference between this 
type of spatial dissipative couphng and bihnear momentum couphng of a system coupled 
to a bath 



have shown that bilinear momentum 
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coupling leads to stochastic acceleration [12, 113| without any violation of the second law of 
thermodynamics |l^. In contrast to spatial coupling which effectively increases the mass 
of the system, momentum coupling reduces it and in the limit of an "Ohmic" coupling, 
the effective mass goes to zero. Hence, the system can undergo stochastic acceleration. 
This observation indicates that perhaps momentum coupling can lead to some rather anti- 
intuitive quantum mechanical results. If it reduces the effective mass, it should amplify 
quantum effects rather than destroy them. This is the topic of this paper. 

In Section II we study the classical and quantum dynamics of a harmonic oscillator 
bilinearly coupled to the momentum of a harmonic bath. We find that as already noted 



in a different context by Cuccoli et al [15|, here too the effective mass is reduced such 
that increasing the momentum coupling increases the thermal variance of the position of 
the quantum particle instead of decreasing it. In Section III we study the dynamics of a 
parabolic barrier. Momentum coupling increases the thermal flux of the particle across the 
barrier as compared to the thermal flux in the absence of coupling. Most interestingly, the 
bath increases the magnitude of the normal mode parabolic barrier frequency, implying that 
it becomes thinner and therefore the tunneling flux through the barrier increases. We show 
that this is indeed the case both above and below the crossover temperature, which now 
increases with increasing coupling strength. We end in Section IV with a Summary. 

II. MOMENTUM COUPLING AND THE HARMONIC OSCILLATOR 
A. Preliminaries 

Our model is that of a harmonic oscillator with mass weighted momentum P, coordinate 
Q and harmonic frequency Q interacting bilinearly with a harmonic oscillator heat bath 



through the momentum ll|. The Hamiltonian then takes the form: 

N N 
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2 2 



(1) 
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where Pj,Xj, j = 1,...,N are the mass weighted momentum and coordinate of the j-th 
bath oscillator whose frequency is uj. The d/s are the bilinear coupling coefficients to the 
particle's momentum. 

Before considering the dynamics of this Hamiltonian, it is appropriate to put it into the 
context of previous studies of dissipative systems. The coupling of the system to the bath 
through the momentum of the bath has been studied previously in a variety of contexts. 
In an early paper, Leggett 16| considered the possibility of two coupling terms, taking 
the form QYl!j=i'^jPj + PJ2f=i^j^j- He then distinguishes between normal dissipation, 
where the Langevin equation of motion is derived for the spatial system coordinate Q and 
anomalous dissipation, where the Langevin equation of motion is derived for the spatial 
system momentum. Neither of these describe the model Hamiltonian given in Eq. [T] above. 
In our model, the coupling term has the form P Ylf=i djPj-, leading to qualitatively different 
dynamics. 

A different model has been considered by Ford et al In their case the coupling to 
the bath takes the quadratic form {^P — Yl!j=i djPj^ ■ It describes the physics of black body 
radiation in which the momentum of the particle is coupled to the magnetic field of the 
vacuum radiation. Such coupling also differs from that given in Eq. [TJ The counter term 
appearing in the model of Ford et al causes a coupling between the bath modes themselves 



and allows by a change of variables 
in form to the standard dissipative 



171 ] to recast the problem into one which is equivalent 



iamiltonian studied in detail in Refs. 
More recently, Cuccoli et al [15|] have studied a momentum coupling model which is 
identical to Eq. [H They term this model as anomalous dissipative coupling. However, as 
already discussed above, the dissipation of this model differs from the one studied by Leggett 
[igI . To distinguish between the two, we have used the terminology momentum dissipation 
for Eq. [TJ 

The formal solution of Hamilton's equations of motion for the j -th bath oscillator is 11| 

xfO) d- /■* 

Xj{t) = Xj{0) cos{ujjt) + sm{uj jt) - — / dt'P{t') sm[ujj{t - t')] . (2) 

UJj UJj Jq 

The equations of motion for the particle are 

N 

Q = P-J2 d,xj (3) 
P = -fi'Q (4) 



Eq. m together with: 



MQ{t) + n^Q = Mfpit) - / dt'P{t')M^p{t - t' 



(5) 



provide a generahzed Langevin equation description for the motion of a particle with the 
effective mass 

N 

M={l + Y.d])-\ (6) 

i=i 

The noise is represented by the momentum fp random acceleration: 



N 



fp{t)=Y,d,u:] 



x(0) 

Xj(0) co?>{ujjt) H — sin(u;jt) 



(7) 



which has zero mean. Its correlation function is: 

N 



P{fp{t)fpm = Y,d]u]cOSUJ,t = 7^p{t) . 

J=l 



(8) 



The brackets denote averaging with respect to the thermal distribution (e ^^). Finally, in 
Eq. we also used the notation: 

N 



ipp{t) = / dt'rjpjt') = d'^-Uj sin (ujjt) 

"^0 7=1 



(9) 



The solution of the generalized Langevin equation ([5]) may be obtained by means of 
Laplace transformation, f{s) = dte~'^^f{t). Using the relation (pp{s) = fip{s)/ s one finds 
from Eq. ([5]) that the Laplace transform of the particle's coordinate is 



Q{s) 



g(o) + sQ{Q) + fp{s) g(o) + sQ(o) + fp{s) 



fi2 



M s 



Noting that 



M (3M 

one readily finds that the classical velocity correlation function is: 



(10) 



(11) 



g(t)Q(o) =i.i.t. -■- 



M s 



M s 



(12) 



where i.l.t. stands for "inverse Laplace transform'' 
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To express results in the continuum limit it is useful to define a momentum spectral 
density as 

N 

Jp{u) = ^-Y,dy,5{u-u,) (13) 
i=i 

where 6{x) is the Dirac function. As a result, the momentum function rip{t) may be 
expressed in terms of the spectral density as: 

r^p(t) = - /" du^^ cos {ut) . (14) 
B. The normal mode transformation 

Additional insight as well as solution of the associated quantum dynamics is facilitated 
by considering the normal modes representation. The Hamiltonian given in Eq. ([1]) has a 
quadratic form and so may be diagonalized. For this purpose we define frequency weighted 
coordinates and momenta as 

Q = VlQ, P = P/n (15) 
Xj = ujjXj, pj=pj/ujj, j = l,...,N (16) 

so that the coordinate part of the Hamiltonian has unit frequency: 

N . , X 2 



2 



/ r! \ . 

.7 = 1 J / , = 1 



Xj 



(17) 



The + 1 normal modes and associated momenta are denoted as UjiPypj = 0,...,N such 
that the normal mode form of the Hamiltonian is: 

*=5E(\X (18) 

j=0 

and the Aj's are the normal mode frequencies. This transformation implies that the vector 
of normal mode momenta Py is an orthogonal transformation of the frequency weighted 
momenta such that 




(19) 

where U is an (A^ + 1) x (A^ + 1) orthogonal transformation matrix. 



Following the same considerations as in the appendix of Ref. 18| one readily finds that 
the normal mode frequencies are the + 1 solutions of the equation: 

fi2 



(20) 



The elements of the transformation matrix are then given by: 



' ' Uko, j = l,...,iV, k = 0,...,N 



I .2 \2 



N 



U 



kO 



-1 



(21) 
(22) 



By considering the 00 element of the (A^ + l)x(A^ + 1) matrix (T" + s^I)~^, (where T" is the 
matrix of second derivatives of the kinetic energy of the Hamiltonian with respect to the 
frequency scaled momenta) one finds the important identity: 



u 



N 

^ S2 + A2 

3=0 J 



1 f}p{s) 



(23) 



This identity then leads directly to all classical results of interest in the continuum limit. 
For this purpose we also define a normal mode momentum function as: 



TV 



j=0 

A spectral density of the normal modes is then defined as [20|: 



(24) 



N 



Tf'(^) = 2 E ^?oA, [S (A - A,) - 5 (A + A,)] 

j=0 



(25) 



One now notes that the Laplace transform of the normal mode momentum function may be 
expressed directly in terms of the original momentum function: 

1 -1 



K s) 



N 2 



• n + \' 



Vp{s) 



M 



(26) 



Using the Fourier decomposition of the Dirac S function (vr5(A) = dt cos{Xt)) we find 
that the spectral density of the normal modes may also be expressed in the continuum limit 
as: 



T(A) = Re AK(iA) 



(27) 
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We further note that at equihbrium: 

(y|> = kBT, J = 0,...,N (28) 

(y]) = \' (< ) = ^'^bT, J = 0, N . (29) 

With these prehminaries it becomes straightforward to solve for thermal correlation func- 
tions. Since the Hamiltonian is diagonal in the normal modes one has trivially that the 
solution of the j-th normal mode is: 

y,{t) = y,{0) cos(A,t) + ^ sin(A,t), j = 0, N. (30) 

The system coordinate is just a linear combination of the normal modes 

Q(t) = = (31) 



SO that 



Similarly 



Q{t)m) = -^- (33) 



These properties of the normal mode transformation become very useful also when 
considering the barrier crossing dynamics, as described in Section llllt below. 



C. Quantum dynamics 

In the quantum regime one may work with the normal modes and use the quantum 
instead of the classical expressions fl28|) . fl29|) . For the equilibrium variances the procedure 
is then as follows. From the quantum mechanical expressions 

(J,|> = ^^coth(5^) (34) 

it follows that 

(<3''> = 4i:"?oA,»th(5fi). (35) 



Using the decomposition [2l| 



w N 1 2a; ^ 1 

COtxifTTX) = 1 > — 



we arrive at 



fc=i 



N oo 

1 + 2 Vm^ 

j=0 k=l J 



(36) 



(37) 



with the Matsubara frequencies Uk = 2'Kk/{h(3). Interchanging the sums, using the fact 



^^qMjo = 1, and the identity (l23l) gives us 



1 + 2E 



02 ( J Vpi>^k) 



(38) 



It is worthwhile to derive a similar expression for the momentum. In this case: 



N 



3=0 



(39) 



with = (n/2Aj)coth(n/3Aj/2). It follows that 



- TV 2 / 

.7=0 ^ 



1 



oo 



j=0 A;=l J 



(40) 



where we used the identity (obtained from Eq. [23] in the limit that s — > 0) that 

1 



N 2 

w|o _ _ 

j=0 3 



E 



(41) 



Finally using again identity fl25]l one obtains 



l + 2fi2^ 



(42) 



It is instructive to derive the quantum mechanical correlations along an alternative route, 
which for spatial friction has been discussed in Ref. QQ a., o., exploits f—al 
principles of quantum statistical mechanics. Namely, since the equations of motion for 
the Heisenberg operators Q{t), P{t) are hnear, the following is true: (i) Due to Ehrenfest's 
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theorem the quantum mechanical averages obey classical equations of motion, (ii) correlation 
functions can be obtained from the quantum version of the fluctuation dissipation theorem, 
and (iii) mean values and second order correlations completely determine the quantum 
dynamics since all random forces are related to a stationary Gaussian process. 

When using (ii) one has to take into account that according to and in contrast to 
spatial friction the time derivative of Q is not identical to P. Hence, to apply the fluctuation 
dissipation theorem we do not start from the equation of motion in position (Q, but from 
the corresponding expression in momentum, i.e., 

im) + ^ {P{t)) - dt'{P{t'))Mt -t') = (43) 

and calculate according to (i) the classical response 

(P(t))= f dt' xp{t-t')F{t'), (44) 

to an external force F{t) applied for t > 0. In Fourier space the above equation reads 

(P(cj)) = Xp{.^) F{uj) so that 

with (pp{ijj) = fjp{—iuj)/{—iuj). 

According to (ii) it is now the symmetrized momentum correlation Sp{t) = 
(l/2)(P(t)P(0) + P(0)P(t)), which is related to the imaginary part of the response function 
Xp = X'p + iXp via 

Sp{uj) = hcoth{ujhl3/2) xp(cu) . (46) 

Further, due to flH|) the anti-symmetrized momentum correlation function Ap{t) = 
(l/'i)[P(t), P(0)] is related to the response function via 

Xp{t) = -le{t) Ap{t) (47) 

with the step function ${■). Thus, in the time domain we arrive at the general expressions 

Sp{t) = — du xp(^) coth(cj;i/5/2) cos(cjt) (48) 
27r i_oo 

and 

Mt) = -7r Xp{^) sin(a;t) , (49) 
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from which following (iii) all real-time correlations can be derived, e.g. {Q{t)Q{0)) = 
(P(t)P(0))/f2^ = [—Sp{t) — iAp{t)]/Q'^. For analytical calculations it is sometimes more 
convenient to work with representations based on Laplace transforms, namely, 

_h 



Ap{s) 
Sp(s) 



Xp[s} 



where we used that Xp{s) 

By considering Sp{0) 
mentum 



Xp{is). 
- lim»_,, 



IXpis) - XpiWnl)] , 



(50) 



sSp{s) one obtains the equilibrium variance in mo- 

2 °o 



+ Z/2 

' n 



(51) 



^ n=-oo ^ ^ n.=l ^ ' V A/ i/„ ' 

which is of course, identical to Eq. fj42l) . For the position = — S'p(0)/r2'^ one must be 

careful when taking the limit in 5'p(0) = lim<j^oo s^S'p(s) due to singularities which must be 
properly subtracted. One gets 



2 ^ 



A/ 



Af 



(52) 



which is identical to Eq. (138|) . When comparing the above expressions with those derived 



for spatial friction 



281], one observes that they can be related to each other by VlQ P, 



which, of course, is a direct consequence of the symmetry of the Hamiltonian ([T]). This means 
though, that all findings known for spatial friction and e.g. ohmic damping can be directly 
translated to momentum friction. We will discuss explicit results in the next Section. 



D. An example 



To obtain a feeling for the various results presented in this Section it is worthwhile to 
consider the specific case of a momentum density with a cutoff: 



TT 



(53) 



where 6{x) is the Heaviside function. One then finds that the Laplace transform of the 
momentum function is 



Vp{s) = 



UJr. 



s arctan 
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and the mass factor 

M = (l + 7cue)-' 

explicitly demonstrating that in the limit that the cutoff frequency goes to infinity, the 
effective mass goes to zero. 

Now, let us first look at the friction dependence of the Laplace transform Xp{^) which 
directly provides the momentum variance and the functions Ap{t) and Sp{t). From 

_ _ Vp{ ) _ _|_ arctan(c<Jc/ s) (54) 
we see that in the limit cUc — > the friction term in xp combines with the potential term 



to yield the effective frequency ^l^s = + 'jcoc of an undamped harmonic oscillator. The 
corresponding variance (P^) as well as Ap{t), Sp{t) are thus trivial. In the opposite limit of 
very large Uc ^ ^,7^^ one has in xp the expression + + f2^77rs/2, which coincides 
with the form of the response function for spatial friction in the ohmic case (friction constant 
7s) with the translation 7<j = r2^77r/2. Thus, the corresponding results can be read off from 
the literature [2], 0]. In particular, for large cutoff and 7i7^7r/2 > 1 the time dependent 
correlations decay to zero monotonously, while for 7n^7r/2 < 1 they decay via damped 
oscillations. In contrast to Ap(t) which displays only classical dynamics, Sp{t) contains an 
additional contribution depending on the Matsubara frequencies that becomes relevant at 
lower temperatures. In the limit of zero temperature Sp{t) then decays in time no longer 
exponentially but algebraically oc Further in the limit of zero temperature one has for 
the momentum variance 



In essence, momentum friction gives for momentum correlations the same results as spatial 
coupling for position correlations. The same is true when the ratios ujc/^l, Uc/'y^l'^ decrease, 
but are still large, and one compares momentum friction with 7, ujc with the well-known 
Drude model for spatial friction, Jp, = 7sC(jexp(— cj/cJc), with 7^ = 7fi^7r/2. 

Let us now turn to (Q^), for which the limit uJc —>■ 00 cannot be taken. We gain 

[i^ )/^L + /5 2^^u^ + n^ + f]2^zy„arctan(a;,/z/J ' ^^^^ 

which for zero temperature, where the sum over Matsubara frequencies must be replaced by 
integrals according to {2/h(3) ^/(t'n) (1/^) / di^fi^)^ reads for strong friction (but still 
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/cb t/ m T/ m 

FIG. 1: Position variances (left, scaled with h/Q) and momentum variances (right, scaled with HQ) 
for momentum friction (solid) and for the corresponding spatial case (dotted) for various values 
of the friction strength 7^1=0.5, 2, 7: left, solid from bottom to top; left dotted from top to 
bottom; right, solid from top to bottom; right, dotted from bottom to top. The cutoff frequency is 
uJc/^ = 10. For the spatial case a model with Drude damping (cutoff lOc) is used with equivalent 
friction strength 7^ = 7il^7r/2. 



(Q^-^ln^^j. (57) 

When dealing with spatial friction, the momentum variance for a Drude model with 7^ = 
'jQ'^tt/2, is given by (P^)spatiai ~ {hQ'^'y /2)ln{uJc/^) [2]. One thus notes that the role of 
the momentum and position variances are interchanged when considering the momentum 
coupling model. In Fig. [1] the position and momentum variances are shown for momentum 
coupling with the spectral density (!53l) together with spatial coupling with a Drude spectral 
density, Jd = 7si^exp(— cj/cUc). Note that the latter one is, apart from its widespread 
use, a sort of minimal model for spatial friction leading to well-behaved variances. In case 
of momentum coupling the spectral density (!53|) plays a similar role. Both models give 
identical results for sufficiently large cutoff and 7^ = fi^77r/2, i.e. f2^(P^)inomcntum(7) = 
(<5^)spatiai(^^77i"/2). They differ though in the low temperature range for the variances 
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momentum > (-P^)spatiai, which are only well-behaved for finite cutoff. Of course, when for 
spatial friction a spectral density is taken that produces the same damping dependence 
in the response function as the one given by we would obtain fully identical results. 
Accordingly, position variances are enhanced for momentum coupling and suppressed for 
spatial coupling and vice versa for the momentum variances. One also observes that at low 
temperatures the momentum variance for momentum friction is suppressed compared to 
the position variance for spatial friction within a Drude model [see f l57|l ]. Accordingly, the 
uncertainty product for T = 

(0^)„(p.)„4,fi-),f^l (58) 



\7ijn^ J 



is smaller compared to the spatial case (Drude damping) by the factor 
[1 - ln(7^7fi/2)/ln(cUe/^^)]. 



III. CLASSICAL AND QUANTUM RATE THEORY 



A. Classical rate theory in the presence of momentum coupling 



We first consider the case of a parabolic barrier Hamiltonian: 

N N 



H 



i=i i=i 



2 2 
UJjXj 



(59) 



In conventional Transition State Theory (TST) one uses the system coordinate as the re- 
action coordinate and the dividing surface is taken to be perpendicular to it. The thermal 
unidirectional flux through the dividing surface is 

Ftst = I dPdQ J] dpjdxje-^''5{Q)Qe{Q) = —j= J] ( ^ ) • (60) 



22 



23|: 



27r 



It is noteworthy that the thermal flux through the dividing surface is larger by the factor 
l/\fM as compared to the conventional TST flux in the presence of only spatial coupling 
to the harmonic bath (Fq = VmFtst)- Momentum couphng causes an enhancement of the 
reaction rate. 

The minimal unidirectional thermal flux is obtained by transforming the Hamiltonian to 
normal modes. Due to the negative force constant associated with the barrier, one will now 



14 



find after diagonalization stable modes witli frequencies Xj and one unstable mode with 
barrier frequency A^. The VTST flux is obtained by considering the flux perpendicular to 
the unstable mode and one finds that: 

N 

I 

FvTST 



By considering the determinant of the second derivative matrix of the kinetic energy in the 
normal mode representation and in the original coordinates one readily finds the identity 

TV N 

At2jjA^2 = fi^2-Q^2_ (g2) 



As a result the ratio of the VTST flux to the TST flux is 

^ = v^|. (63) 



The analog of the Kramers-Grote-Hynes equation |24l . l25l ] for the normal mode barrier 
frequency in the presence of spatial diffusion is obtained from Eq. ( l20l) . except that one 
substitutes the stable mode frequency with the unstable mode frequency. After a bit of 
rearranging one finds: 

(a4) 

from which it follows that -\/M ^ < 1 and as expected the VTST flux is lower than the TST 
flux. Note however that Eq. (!64l) may also be rewritten as: 

showing that momentum coupling leads to a thinning of the barrier, instead of the usual 
broadening of the barrier found as a result of spatial coupling. In fact the ratio ^^^^"^ = 
implying that even the minimal VTST flux is larger then the flux in the absence of 
coupling. This is proof that classically, momentum coupling leads to an increase of the 
thermal parabolic barrier crossing rate, as compared with the absence of coupling. 

It is instructive to consider a specific example, namely the spectral density given in 
Eq. ( l53l) . Using the reduced values x = X^/Q^,Wc = cUc/fi"'" and g = we plot in Fig. [2] 
the reduced barrier height as a function of the reduced momentum coupling coefficient g for 
three representative values of the reduced cutoff frequency Wc- As noted from the figure, 
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FIG. 2: Classical enhancement of the reactive flux as a result of momentum coupling to a harmonic 
bath: The parabolic barrier frequency is plotted as a function of the coupling strength for three 
different values of the cutoff frequency. The reduced barrier frequency is also the ratio of the exact 
thermal unidirectional flux to the flux in the absence of coupling to the bath. The solid line is for 
Wc = 10, the dashed line for Wc = ^ and the dotted line is for Wc = 0.5. 

only when the cutoff frequency is large does one get an appreciable increase in the barrier 
frequency. In the limit oi g ^ Wc one has that x ~ ^JgWc- In Fig. [3] we then plot the 
ratio FvTST / Ftst for the same parameter range as in Fig. O One notes that there is an 
appreciable effect on the transmission factor only when the cutoff frequency is much larger 
than unity. In contrast to the spatial coupling case, here the transmission factor tends to 
unity both in the weak and strong momentum coupling limits. 

B. Quantum rate theory in the presence of momentum coupling 

Given the fact that momentum coupling causes a decrease in the barrier width, it is 
interesting to study whether it could then cause an increase in the tunneling probability, 
since the tunneling probability is exponentially sensitive to the width of the barrier. We will 
consider two cases. The first is transmission through a parabolic barrier. The second will be 
consideration of tunneling through an anharmonic barrier in the limit of weak momentum 
coupling and at temperatures which are below the crossover temperature between tunneling 



16 



1.0 




Q 5 I I I I I I I I 

■ 5 10 15 20 

9 

FIG. 3: VTST solution for the classical transmission factor through a parabolic barrier potential as 
compared to choosing the system coordinate as the reaction coordinate. The three lines correspond 
to the same values of the cutoff frequency as in Fig. [2l 



and thermal activation. We will see that the two limits give qualitatively identical results. 
Above the crossover temperature, where the parabolic barrier approximation is valid we 
find that momentum coupling indeed causes an increase of the rate which is even greater 
than the classical enhancement. This is due to the thinning of the parabolic barrier. At 
low temperatures, below the crossover temperature, momentum coupling leads even to an 
exponential increase in the tunneling rate. 



1. Above the crossover temperature 

In the limit of jx parabolic barrier, we follow the same reasoning as for spatial diffusion, 
as given in Ref. 
the barrier is 



a]. In the absence of momentum coupling, the thermal fraction through 



while the VTST tunneling fraction is given by the same expression, except one must replace 
the bare frequencies everywhere with the normal mode frequencies. One thus has that the 
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quantum transmission factor is 

However we know that the barrier frequency > Q^. The denominator with the sin 
function will diverge at the temperature h(3c\^ = 2tt, that is the crossover temperature 
will be greater than the crossover temperature in the absence of coupling. In contrast to 
spatial coupling, which reduces the crossover temperature, momentum coupling increases it. 
At the crossover temperature, the quantum transmission factor will diverge, implying that 
the quantum transmission factor is indeed larger than the classical one. As noted, for the 
parabolic barrier, momentum coupling enhances tunneling. 

The result for the tunneling fraction has to be transformed so that it may be expressed 

in the continuum limit. For this purpose one uses the infinite product representation of the 

! 

sin and sinh functions (as detailed in Ref. [2fi]) as well as the identity: 

det (T" + .^I) = (-A^^ + s^) n (A? + s^) = ( + s^ + 1^^^^ U + s^) 

j=i ^ ^ j=i 

where as before T" represents the matrix of second derivatives of the kinetic energy with 
respect to the (frequency scaled) momenta. The last equality on the r.h.s is obtained by car- 
rying out explicitly the evaluation of the determinant, using the frequency scaled momenta, 
as given in (ITS!) and (fT6|) . One then finally finds that 

11 ,,2 f^t2 , nt2'7p(^'fc) ~ OX^P ^^^^ 



where = -j-^ are the Matsubara frequencies and Sp is the Wolynes factor ^] associated 
with momentum coupling. 

For the specific spectral density given in (!53|) one finds that the Wolynes factor is: 

= n — — (^0) 

ti yl - (l + l^k arctan (^^ j j 

showing clearly that the quantum enhancement of the rate due to momentum coupling is 
greater than the classical enhancement, since in the classical limit the Wolynes factor tends 
to unity. In Fig. H] we plot this Wolynes factor as a function of the reduced momentum 
coupling coefficient g = jfl and reduced cutoff frequency Wc = uid^- 
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FIG. 4: Macroscopic enhancement of the quantum tunneling rate above crossover due to momen- 
tum coupling to a harmonic bath. Panels (a)-(c) correspond to the cutoff frequencies 10,4 and 
0.5, respectively. In each panel we plot the dependence of the Wolynes enhancement factor as a 
function of the coupling strength for three different temperature values. Defining 6 = hf3Q^, the 
dotted, dashed and solid lines correspond to ^ = 0.5, 1, 3 respectively. Note the change in scale 
for the coupling strength as one goes from panel (a) to panel (c). For low cutoff frequencies, the 
enhancement is weakened. 
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2. Tunneling through an anharmonic barrier below the crossover temperature 



For tunneling through an anharmonic barrier potential V{Q) with a barrier top located 
at Q = Qi, the escape rate F is most conveniently calculated from the imaginary part of 
the free energy and thus from the imaginary part of the partition function of the unstable 



system 



0,3- 



In the path integral representation one has for the reduced system 

Z = Tr{e-^^} = ^V[Q] e'^^^ 1^1/'^ , (71) 

where one sums over all periodic paths running in the imaginary time interval hjS through 
the inverted barrier potential —V{q). Here for momentum friction the effective action is 
found to read 

SM = J^ '^^[2^^''^' + ^^^^ + 2 J, daQ{T)M{r-a)Q{a)\ (72) 



with the kernel given by 



n 
2 

71 



For high temperatures when z/„ — 00 one has M(r) =: 5{t) : M^jd'j with M = 
1/(1 + '^jd'j). Combining the friction term with the bare kinetic term then leads to an 
effective kinetic term of the form so that friction appears simply as an effective 

mass. Accordingly, the thermal activation factor F oc exp[—(3V{Qb)] is independent of fric- 
tion since it is determined by the constant path at the barrier top Q = Qh, while dissipation 
influences the prefactor as specified above. The same is true for lower temperatures above 
the crossover when quantum fluctuations in the rate prefactor lead to an even stronger 
increase of the escape rate as shown in the previous Section. 

Below the crossover temperature the contribution of the bounce orbit Qb dominates the 
imaginary part of the partition function and the exponential factor in the rate contains its 
action Sb = S[Qb], i-e. F oc exp(— Ss//?,). Now, let us consider weak friction. In this case 
we write for the bounce path Qb = Qo + ^Q, where Qo is the bounce path in absence of 
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dissipation, which obeys Qo — V'{Qo) = 0. Upon inserting Qb into the effective action (172 
one finds that to lowest order in the friction one has Sb = S[Qo] + A^o with 



1 fH/B fhfi 

ASo = - dr daQoiT)M{T-a)Qoia). (74) 
^ Jo Jo 



This correction can be easily expressed in Fourier space by using Qq 
{l/h(3) XlnQn''' 6xp(zz/„r), wherc one may choose the phase of the bounce such that Qn^ 
Q-n = Qn'' ■ This way, we find with m„ = m_„ and Un = —i^~n th.at 



AS, 



;^j:\Q'n'\'-l^n. (75) 



2hp 

Apparently, A^o < so that Sb < Sq meaning that the probability for quantum tunneling 
is exponentially enhanced due to momentum friction. Physically, since the effective mass of 
the combined kinetic terms M(r)+ : 6{t) : has Fourier components that are always smaller 
than 1, the particle's kinetic energy becomes smaller relative to its potential energy. Thus, 
for a dynamical orbit like the bounce path the action decreases. In contrast, spatial friction 
leads always to a rate suppression since effectively it provides an additional contribution to 
the potential energy. For stronger coupling to the heat bath the bounce orbit can easily be 
calculated numerically along the lines described in 28 1. 



IV. DISCUSSION 



In this paper we analyzed the influence of a harmonic heat bath coupled to a system 
with potential energy via a bilinear interaction between the system's momentum and the 
individual momenta of the bath degrees of freedom. For a harmonic oscillator and momen- 
tun, d.sipat„n we End . *o noted ,ua„tat,vely by CnccoU et a, Q that the coapHn. to 
the bath increases the delocalization of the position of the oscillator. The thermal variance 
increases with increasing coupling, exactly the opposite of the behavior found for spatial 
coupling. 

For escape over a barrier, in the high temperature regime, where thermal activation pre- 
vails, momentum friction leads to an increase of the flux across the barrier. At somewhat 
lower temperatures quantum fluctuations come into play and enhance the flux even more. 
Above the crossover temperature between tunneling and thermal activation, friction appears 
only in the prefactor of the rate expression. Below the crossover temperature, where quan- 
tum tunneling dominates, it reduces the action of the bounce path and thus exponentially 
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increases the decay rate. In contrast to the case of spatial friction, the crossover tempera- 
ture increases with increasing friction so that for pure momentum couphng quantum effects 
would be observable even in the high temperature domain. 

In this paper we considered the case of pure momentum coupling. Any realistic system 
will be influenced by both, momentum and spatial coupling. The crucial issue is then to what 
extent the latter one suppresses the impact of the former one. Specific systems to be studied 
in the future in this respect include molecular compounds, mesoscopic islands coupled to 
fluctuating charges and the transport of charged particles moving under the influence of 
random magnetic fields [12I, Q • 
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